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Abstract 

In the Large-Nc limit of QCD, two-point functions of local operators become Har- 
monic Sums. I review some properties which follow from this fact and which are relevant 
for phenomenological applications. This has led us to consider a class of Analytic Number 
Theory Functions as toy models of Large-Nc QCD which I also discuss. 
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I Introduction 

Many of us would like to know the answer to the following questions: 

• What is the effective Field Theory of QCD at Long Distances ? 

• How does QCD fix the couplings of the Chiral Lagrangian of the Nambu-Goldstone 
modes of the spontaneously broken chiral-SU(3) flavour symmetry ? 

• Can we answer these questions, perhaps more easily, within the framework of QCD in 
the Large-Nc limit [T] ? 

In that respect, it has been shown H that if the confinement property of QCD persists 
in this limit there is also spontaneous chiral symmetry breaking and the Hadronic 
Spectrum consists then of an infinite number of Narrow States [7]. 

Unfortunately, in spite of the successes of the Standard Model, the answer to these questions 
remains unknown. What I shall do here is to provide a few comments related to them. 

I.l General Comments 

1. Independently of the Large-Nc approximation, the couplings of the effective Chiral 
Lagrangian of the Strong Interactions of the Nambu-Goldstone modes, can be identified 
with the coefficients of the Taylor expansion of appropriate QCD Green's Functions. 

2. By contrast, most of the couplings of the effective chiral Lagrangian of the Electroweak 
Interactions of the same Nambu-Goldstone modes are given by integrals over all the 
range of euclidean momenta of appropriate two-point functions with soft insertions of 
local operators. Their determination, therefore, requires a precise matching of the short 
distance and the long distance contributions to the underlying QCD Green's functions. 

Typical terms of the chiral Lagrangian are: 



£eff = tr (i?^[/D^[/t) +Cwzw + Lio tr (u^ Fr^.U F^"^ +■■■ 

TTTT— )-7r7r , K^Tvey ir—^evy 

e^C triQnUQLU^) - %KdK*s 98 (d^UD^uA +■■■ . (1.1) 

_e2C-^(7r+7r-+ K+R-) ^ ' 

f4 K^tttt, K^tt-ktt 



Here U denotes the 3x3 unitary matrix in the u,d,s flavour space which collects the Nambu- 
Goldstone fields and which under chiral rotations transforms as C/ — )■ VrUV^; D^U denotes 
the covariant derivative in the presence of external vector and axial-vector sources. The first 
term in the first line is the lowest order effective Lagrangian in the sector of the strong inter- 
actions 15 with the pion-decay coupling constant in the chiral limit where the light quark 
masses u,d,s are neglected (-F^r — 90 MeV). The second term stands for the anomalous Wess- 
Zumino-Witten [9l[T0] effective Lagrangian of 0{p'^). The third one shows a typical term of 
0{j)^) in the chiral counting |llj . with Lio a coupling constant which is not fixed by symmetry 

^See refs. [H |31 IH [S] . For a lucid exposition see ref. [B]. 
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requirements alone. The first term in the second hne corresponds to the lowest order effective 
Lagrangian, which is [©(p*^)] in the chiral counting [12], and which appears when photons are 
integrated out in the presence of the strong interactions {Ql = Qr = diag.[2/3, —1/3,-1/3] 
and e is the electric charge). The second term in the second line is the lowest order effective 
Lagrangian in the electroweak sector which induces non-leptonic X-decays @. Typical physi- 
cal processes to which each term contributes are indicated under the braces. Each term in the 
chiral Lagrangian is modulated by a coupling constant (F^ , Lio , • • • C , ■ ■ ■ in Eq. (jl.ip ). 
which encodes the underlying dynamics responsible for the appearance of the corresponding 
effective term. The evaluation of these couplings from the underlying QCD Lagrangian is the 
question we are interested in. 

1.2 The Left Right Correlation Function as an Example 

As a precise example of the comments in the previous subsection let us consider the left -right 
correlation function: 

n{^^(g) = 2i j (i^2;e^«-^(0 I T (L'^(x)i?^(0)t) | 0) , (1.2) 
with left and right currents: 

L^(x) = fi(x)7^^(l - -iMx) and ii^(x) = (i(x)7^^(l + 75)w(x) . (1.3) 

The discussion here, unless explicitly mentioned, does not use the Large-Nc approximation. 
In the chiral limit where the light quark masses are set to zero, this two-point function only 
depends on one invariant function (Q^ = —q^ > for q'^ spacelike) 

KUq) = ii'<f - b'^'q') nLR(Q') . (1.4) 

The self-energy function HhRiQ^) in the chiral limit vanishes order by order in QCD pertur- 
bation theory and is an order parameter of the spontaneous breakdown of chiral symmetry 
for all values of the momentum transfer [T3|. In what follows we shall be working in this 
limit. 

The Taylor expansion of ^lr{Q^) at low is a power series in Q^: 

- Q'nLR(Q2) = f2 + 4Lf,Q^ + OiQ^) , (1.5) 

More precisely 

L\q = Lio(^) + Goldstone one loop corrections (1.6) 

where Lio(/i) denotes the 0{p'^) coupling in Eq. (jl.ip renormalized at the scale fi. In the 
1/Nc-expansion, the Goldstone loop corrections are subleading and, therefore, only the tree 
level couplings survive in the Large-Nc limit, which is one of the important simplifications of 
this limit. 

The function nLR(Q^) also governs the coupling constant C in Eq. (II. ID and gives a mass 
of electromagnetic origin [15j (and also from the integration of the heavy Z electroweak 
boson [16]) to the vr^ and particles: 

2e^C n n a 



2 2 

m_± = m 



K^ = ^^fdQ' [-Q'U^Q')] . (1.7) 



^Gf is the Fermi constant and Vud, V^g matrix elements of the flavour mixing matrix. The couphng gs_ 
governs the strength of the dominant A/ = 1/2 transitions. For more details see e.g. ref. [13] and references 
therein. 
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This integral converges in the ultraviolet because [T7] 

lim HMQ') - O (M!) . (1.8) 

Furthermore, it has also been shown |18l [T9] that 

- Q^nLR(Q^) > for aU < Q2 < oo , (1.9) 

which in particular ensures the positivity of the integral in Eq. (jl.7p and thus the stability of 
the QCD vacuum with respect to small perturbations induced by electroweak interactions. 

II Large Nc Models and Phenomenology 

QCD models of spontaneously chiral symmetry breaking, like the Constituent Chiral Quark 
Model (CxQM) of Georgi and Manohar [201 [2T1 [22l [23] , or the more sophisticated Extended 
Nambu-Jona-Lasinio model (ENJL) [241 [25] have been rather successful in reproducing the 
phenomenological determinations of the couplings of the chiral Lagrangian in the strong 
interaction sector. However, they fail in general to provide the required matching between 
short and long distances which is needed in order to evaluate the couplings of the chiral 
Lagrangian induced by the Electroweak Interactions, i.e. couplings like C and in Eq. (jl.ip . 
Because of this, in the phenomenological applications, these models have progressively been 
replaced by a more direct approach where the relevant Green's functions are approximated 
by a finite number of the Large-Nc QCD hadronic spectrum of narrow states. In fact, the 
Resonance Chiral Lagrangians of the type discussed in refs. |12[ [26] and their extensions (see 
e.g. ref. |27] and references therein), can be viewed as simplified versions of the Large-Nc 
QCD Hadronic Lagrangian when limited to a finite number of states. 

The methodology which has been suggested (see e.g. ref. [28| for a review), consists 
in fixing the couplings and masses of a Minimal Hadronic Ansatz (MHA) of narrow states 
which contribute to a specific Green's function in such a way that, on the one hand the 
short distance behaviour predicted by the operator product expansion (OPE) [17] of the 
underlying Green's function in Large-Nc QCD is satisfied and, on the other hand, the long 
distance behaviour constraints governed by the effective chiral Lagrangian in the sector of 
the Strong Interactions alone are satisfied as well. As an example of this MHA approach, let 
us consider again the integral in Eq. (II. 7|) which fixes the coupling C. The MHA in this 
case requires the presence of three states: the massless pion pole, a vector state with mass 
My and an axial-vector state with mass Ma- The constraint that ^\m{Q^) satisfies the OPE 
at short distances (see Eq. (jl.Sp ) implies that 0{1/Q'^) terms and 0{1/Q^) terms must be 
absent. This, as we shall later discuss, also implies the validity of so called Weinberg sum 
rules, and fixes the couplings of the V and A states in the spectral function of the left-right 
correlation function with the result 

Inserting this function in Eq. (II. 7p gives a prediction for the tt"*" — tt*^ = AttItt mass difference E|: 

Am^ = (4.9 ± 0.4) MeV , (2.2) 

'^This is the result for i^^ = (87 ± 3.5) MeV, My = (748 ± 29) MeV and = ^ = 0.50 ± 0.06. These values 
follow from an overall fit to predictions of the low energy constants in the MHA to Large-Nc QCD. 
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to be compared with the experimental value: Atti-tt = (4.5936 it 0.0005) MeV. 

The MHA approach to Large-Nc QCD has led to a remarkable set of interesting pre- 
dictions for some of the couplings of the Electroweak Lagrangian in the chiral limit (two 
representative references are |29l I13j ). The incorporation of chiral corrections, however, be- 
comes technically rather cumbersome and, above all, the question of the reliability of the 
approximation with a finite number of narrow states to Large-Nc QCD, to which I shall later 
come back, remains open (see e.g. ref. [30j). 

II. 1 Comment on Minkowsky versus Euclidean 

The spectral function of the left-right correlation function defined in Eq. (jl.2p can be obtained 
from measurements of the hadronic r-decay spectrum (vector-like decays minus axial-vector 
like decays). Figure 1 below shows the experimental determination of ;^ImnLR(t), obtained 
from the ALEPH collaboration [31j data at LEP, versus the invariant hadronic mass squared 
t in the accessible region 0<t<ml. This plot is a good experimental proof of spontaneously 
chiral symmetry breaking in Nature. If the symmetry was realized a la Wigner-Weyl, the 
shape should be a straight horizontal line all the way down to zero. Also shown in the 
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Fig. 1 The Spectral Function ilmnLR(i) compared to the MHA to Large-Nc QCD. 

same Fig. 1 is the simple spectrum of the MHA approximation to Large-Nc QCD which, as 
discussed before, consists of the pion pole (not shown in the figure), a vector narrow state 
(the first positive vertical line) and an axial-vector narrow state the second negative vertical 
line) . The heights of these vertical lines represent the strengths of the couplings of the narrow 
states in the spectral function. At this level of approximation, and in the chiral limit, the 
rest of the vector and axial-vector states are degenerate and, therefore, they cancel in the 
difference, a fact which in any case is predicted by QCD at high energies when perturbation 
theory applies, and is reproduced by the asymptotic horizontal continuum line shown in the 
figure. 

Looking at the plot in Fig. 1, one can hardly claim that the MHA approximation in the 
Minkowsky region reproduces the details of the experimental data. This is to be contrasted, 
however, with what happens in the Euclidean region. With HhRiQ^) determined from the 
spectral function by the unsubtracted dispersion relation 
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Fig. 2 Plot of ^^nLR(Q^) in the Euclidean region. The solid curve is the one corresponding to 
the MHA to Large~Nc QCD and the dotted curve the one from the experimental data in Fig 1. The 
other curves are the predictions of the Constituent chiral Quark Model(CxQM) and the Extended 
Nambu-Jona-Lasinio Model (ENJL). 

nLR(Q^)= / dt——.-lmULR{t), (2.3) 
the corresponding plot of tlie function ^-nLR(Q^) versus the Euchdean variable nor- 

^ TV 

malized to vector mass My, as shown in Fig. 2 (the solid curve), reproduces rather well the 
dotted curve, which is the one resulting from the experimental data in Fig. 1. This is in fact 
a generic feature: Green's functions in the Minkowski region have a lot of structure, while the 
corresponding shape in the Euclidean region shows a smooth behaviour. As shown in this ex- 
ample, the simple MHA in the euclidean region already provides a rather good interpolation 
between the asymptotic regimes where, by construction, it has been constrained to satisfy 
the lowest order chiral behaviour and the leading and next to leading OPE constraints. This 
good interpolation is the reason why the integral in Eq. (II. 7p . evaluated in the MHA, already 
reproduces the experimental result rather well. One can also see in Fig. 1 that models which 
fail to incorporate the matching to short distances, like the CxQM and the ENJL model, fail 
to reproduce the shape of the experimental curve already at rather low values of Q^. 

II.2 More on the CxQM 

Concerning the failure of the matching between long and short distances in Models of Large- 
Nc QCD, I would like to point out that there is, however, a class of low-energy observables, 
governed by integrals of specific QCD Green's functions, for which the CxQM predictions 
encoded in the MG-Lagrangian, in spite of its limitations, can be rather reliable. This is the 
case when the leading short distance behaviour of the underlying Green's functions of a given 
observable is governed by perturbative QCD. Interesting examples of this class of observables 
are the decay rate vr'' — t- e^e~, the Hadronic Vacuum Polarization and the Hadronic Light- 
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by-Light Scattering contributions to a low-energy observable like the anomalous magnetic 
moment of the muon: ^{g^ — 2). Furthermore, as recently pointed out by Weinberg [22j . 
the MG-Lagrangian in the Large-Nc limit, modulo the addition of a finite number of local 
counterterms |23j . is a renormalizable Lagrangian. Calculations with the MG-Lagrangian, 
compared to those with the more sophisticated approaches described above, have the advan- 
tage of simplicity and, when applied to this class of low-energy observables, can provide a 
check to the more elaborated phenomenological approaches. The fact that one is dealing with 
a renormalizable Quantum Field Theory is, of course, a welcome feature. 
The CxQM effective Lagrangian in question is the following: 



CcxQuix) = iQr {d^ + + iG^) Q - '-QA Ql^l^i^^Q - MqQQ -^Q (S - 75A ) Q 

V ' 

M-G 



1 



a=l 



M-G 



M-G 



L5 tvDM^D^'U{x^U + U^x) + ^8 HUxWx] + U^xU^x) ■ 



(2.4) 



The underbraced terms are those of the MG-Lagrangian, but in the presence of external 
SU{3) vector and axial-vector a^{x) sources. The field matrix U{x) is the same 3x3 

unitary matrix in flavour space which collects the Goldstone fields as in Eq. (jl.ip . The vector 
field matrix D^U is also the same covariant derivative of U as in Eq. (II. ip with respect to 
the same external sources: 



and, with U = 



1 

2 L' 



(2.5) 



(2.6) 



The gluon field matrix in the fundamental representation of color SU (3) is G^{x) and G^^} (x) 
its corresponding gluon field strength tensor. The presence of external scalar s{x) and pseu- 
doscalar p{x) sources induces the extra terms proportional to 



X = 2B[s{x) +ip{x)] , 



(2.7) 



where B, like is an order parameter which has to be fixed from experiment. When these 
sources are frozen to the up, down, and strange light quark masses of the QCD Lagrangian, 



and then 



X = 2BAi , with Ai = diag(mu , , rus) 



(2.8) 



(2.9) 



With the axial coupling fixed io qa = 1, the extra couplings L5 and Lg are the only 
terms which are needed to absorb the ultraviolet (UV) divergences when the constituent 
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quark fields Q{x) are integrated out |j. If one wants to consider the case where photons are 
also integrated out then, to leading order in the chiral expansion and in the electric charge 
coupling e, the last term in the second line is also required to absorb further UV~divergences. 
Loops involving pion fields are subleading in the 1/Nc-expansion and hence, following the 
observation of Weinberg in ref . [22] , the Lagrangian in Eq. ()2.4p , when considered within the 
framework of the large-Nc limit, is a renormalizable Lagrangian. 



Ill General Properties of Large Nc QCD 

An interesting feature of Large-Nc QCD which I shall discuss next is the fact that two-point 
functions of color singlet composite local operators become simple Harmonic Sums in this 
limit. 

In full generality an Harmonic Sum is characterized by a Base Function (in our case 
the kernel in the dispersion relation which the two-point function in question obeys) and a 
Dirichlet series. 

oo 

E(s) = ^A„Mn% (3.1) 

n=l 

where the are called the Frequencies (the position of the narrow states in the Minkowski 
region in our case) and the the Amplitudes (the residues of the corresponding poles). 
The simplest example of a Dirichlet Series is the Riemann zeta function: 

oo 



as) = Y,n-'= n T^' I^«(^)>1' (3-2) 



n=l primos(p) 

where A„ = 1 and fin = n. The Euler product expression in the r.h.s. extends to all prime 
numbers p. 

III.l The IIlrXQ^) Self Energy as an Harmonic Sum 

Let us consider again the left -right correlation function in the chiral limit as an example to 
illustrate these properties. In the Large-Nc limit, the spectral function ilmnLR(t) is the 
sum of an infinite number of narrow states 

^ oo 

-lmnLR(t) = -F^5{t) + ^(-1)"+1f2 6{t - M^) . (3.3) 

n=l 

with positive weights for the vector-like components and negative weights for the axial- 
vector-like components. The first term is the contribution from the pion pole. Inserting this 
spectral function in the dispersion relation in Eq. (12. 3p we find 



|,n,R(Q^) = l-5:(-l)-^||-i^. (3.4) 



This is a typical Harmonic Sum which we can write as follows: 



q2 oo / j^j2\ 

1 + ^nLR(Q2) = ^ A„ i?LR y iln-^ j , (3.5) 



®We disregard divergent couplings involving external fields alone to lowest order in the chiral expansion. 
^For a clear introduction to the appropriate mathematics literature see e.g. ref. [32| . 
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where (Mp = Mi 



p2 i\j2 

A„ = (-ir+i^ and /xn = v^ (3.6) 



and where the Base Function is 

^lr(^) = ■ (3.7) 

The crucial property of Harmonic Sums is that they have a factorizable Melhn-Transform[f|. 
In our case 



M 

with 



1 + |JnLR(Q2) 



(s) =A^[i?LR(x)](s) S(s), (3.8) 



M[Bi^K{x)]is) = r{s)r{l-s), (3.9) 



and S(s) the Dirichlet series in Eq. 13.11 with A„ and /i„ defined in Eqs. 13.61 Inverting the 
Mellin transform in Eq. ()3.8p results then in the following Mellin-Barnes representation: 

c+joo / _ ■ 



i?2 

TV 



1 /■ / \ 

nLR(Q2) = i-— y ^MgrJ s(.)r(.)r(i-.), (3.10) 



where the integration path lies in the so called Fundamental Strip [32] which is defined by 
the intersection of the convergence domain of the base function (Re(s) G ]0, +1[ in our case), 
with the domain of absolute convergence of the Dirichlet Series S(s). Notice that in full 
generality, independently of the Large-Nc approximation, ^hK{Q^) also has a Mellin-Barnes 
representation like the one in Eq. ()3.10p . except that then S(s) is the Mellin transform of the 
spectral function ilmnLR(t) where the tilde means that the contribution from the pion pole 
has been removed: 



, ^ M^^ dt f t I ~ , , 



This integral, in the Large~Nc limit, becomes the Dirichlet Series: 



The Mellin-Barnes representation is very useful. All the dynamics is encoded in the 
factorized Dirichlet series S(s) or, in full generality, in the Mellin transform of the Spectral 
Function. The dependence in is completely factorized from the details of the spectrum, 
and one can read off the asymptotic behaviours for small (i.e. the chiral expansion) and 
large-Q^ (i.e. the short distance expansion or OPE expansion) in a straightforward way 
These asymptotic behaviours can be obtained from the so called inverse mapping theorem [32] 
as follows: 



®The Mellin transform of a function F(z): M[F(z)](s), is defined by the integral M[F{z)]{s) = 
Jq°° dzz"'^ F{z), in the domain of s where the integral exists and elsewhere by its analytic continuation. 

'■'Applications of this property of the Mellin-Barnes representation in QED and QCD have been recently 
discussed in refs. 1331 1341 1351 and references therein. 
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• The coefficients Sp^^ of the singular series expansions at the left of the fundamental strip 
i.e. 

7V(p) 

r(.)r(i - .) X 5] ^ j-^^ (3.13) 

5Rp>0 fc=0 ^ 

govern the asymptotic expansion: 

• The coefficients bp ^ of the singular series expansions at the right of the fundamental 
strip i.e. 

A^(p) b 

s(.)r(.)r(i-.)x Etjz^ (3.i5) 

SRp>l fc=o ^ 

govern the asymptotic expansion: 

The sum structure over the poles p in the r.h.s. of these equations is completely general 
for any two-point function. When applied to the left-right correlation function there are, 
however, further specific restrictions to take into account: 

• The fact that in the OPE of TlhKiQ^) there are no operators of dimension two and four 
implies that 

ao,o = l, and ai,o = 0, (3-17) 

and hence 

S(0) = 1, and S(-1)=0. (3.18) 

This is precisely the content of the Weinberg Sum Rules, conventionally written as 
follows: 

f°° dt 1 

/ ^-lmnLR(t) = 1 , 1st Weinberg Sum Rule, (3.19) 
Jo F^T^ 

f°° dt t 1 

/ — J— -r-ImnLR(t) = 0, 2nd Weinberg Sum Rule . (3.20) 
Jo Fi 

Notice, however, that in full generality sum rules like 

at the left of the fundamental strip where s < 0, have to be understood as analytic 
continuations of the Mellin transform of the spectral function (or the Dirichlet series 
in the Large-Nc approximation) in the region of s of absolute convergence. It is in this 
sense that the two Weinberg sum rules should be understood. 

Because of the absolute convergence of S(s) at the right of the fundamental strip, all 
the poles in the corresponding series expansion must be simple poles, hence A^(p) = 
in Eq. ()3.16p . which means that the expansion at small is a pure power series. 



and 
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IV Toy Models of Large Nc QCD 



As we have seen, the phenomenological apphcations of Large-Nc QCD we have discussed in 
Section II rely on the hmitation of the hadronic spectrum to a finite number of narrow states. 
In order to test this approximation and in the absence of a solution of QCD in the Large-Nc 
limit, one has to resort to models to investigate the issue 0. As we shall see, some of these 
models bring in functions of Analytic Number Theory with quite interesting features. 

IV.l The Hurwitz Model HlrXQ^) 

This is a model where the spectral function of the left-right correlation function consists of 
an infinite number of narrow states 

^ oo oo 

-lmnLR(t) = -F^6{t) + F2 V 6{t - M2 - na^) - V S{t -M\- na^) , (4.1) 

TT ^ — ' ^ — ' 

n=0 n=0 

with the vector states and the axial-vector states equally spaced by an amount a. Implement- 
ing the two Weinberg sum rules in this spectral function reduces the number of independent 
parameters to F-j^ , the lowest vector mass My and the value of qa = -ttt , with Ma the lowest 
axial-vector states. The other parameters are then fixed as follows: 

F^ = Fl = F^^\±^, a^ = M^(l + l-), and 9A = ^. (4.2) 

Notice that this forces the inequality 

gA<l and therefore M\> My . (4.3) 
The Mellin-Barnes Representation in this model is then: 



c+ioo -, 



Fq 1- qa 2iri J V 9 A 



C — lOO 



c ( s, -c(s, ](^] nsm - s) , (4.4) 



i + 5Ay V ^ + 9aJ j \ Q 



where C{s, v) is the Hurwitz function, a generalization of the Riemann zeta function, defined 
by the series: 

oo ^ 

C(s, v) =y — , Res > 1 , with v a fixed real number < f < 1 , (4.5) 

and its analytic continuation. For w = 1 it reduces to the Riemann zeta function. 

The Hurwitz function is a special case of the so called Dirichlet L- functions [39]. It has 
an integral representation in terms of a Mellin transform: 

/■oo p—tix 

r(s)C(s,v)= dx x'-^ , Res>0, (4.6) 

Jo 1-e ^ 
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Fig. 3 Plot of the Mellin Transform of the Hurwitz-Model in Eq. for qa = 1/2. 



which sets the basis for its analytic continuation [39j. The so defined C{s,v) is analytic for 
all s except for a simple pole at s = 1 with residue 1. 

Comparing Eq. (|4.4p with Eq. (|3.1Up . we see that the relevant Mellin transform of the 
Hurwitz-Model is 



Shm(s) = -r^^— ( 1 + — 
1 - 5A V 9 A 



l-s 



'^ + 9aJ V ^+9A 



(4.7) 



A plot of this Mellin transform is shown in Figs. 3. Notice the zeros at s = — 1, —3, • • • . These 
zeros are in fact quite generic of the model and they occur at odd negative values of s. The 
origin of it is due to the property that, for s = —m with m = 0, 1, 2, . . . : 

A/ N Bm+i{v) 

Q[-m,v) = -— , (4.8) 

m+1 

where Bm+i{v) denotes the Bernoulli polynomial of degree m+1. This fact plus the symmetry 
property 

Bm{l-x) = {-irB^{x), (4.9) 

is at the origin of the zeros of the Mellin transform at s = —1, —3, —5, .... Notice also that 
the values of the Mellin transform at negative even integer values s = —2, —4, —6, . . . are 
fixed by the values of Bernoulli polynomials (m = 0, 1, 2, 3, . . . ): 

.„„,_2„.) . (l + l±(-f B.„,, (-^) . ,4.10) 

Myl-9A\ 9aJ m + 1 \l + gAj 

At this stage, it is interesting to compare the Mellin transform of the Hurwitz model with 
the one of the MHA-approximation discussed in Section II. This is shown in Fig. 4. What 
emerges from this comparison is the fact that the two Mellin transforms do not differ much 



'Early work on models of this type can be found in refs. [3S1 133 EB] and references therein. 
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Fig. 4 Mellin Transforms of the Hurwitz Model (solid curve) and MHA (dashed curve). 



in tiie region —1.5 ^ s < 1.5, the reason being that they both satisfy the two Weinberg sum 
rules. The result that the pion mass difference was rather well reproduced with the MHA 
approximation in Eq. (j2.2p is due to the fact that this observable is governed by the slope of 
the Mellin transform at s = —1: 

The two Mellin transforms in Fig. 4 have, however, very different behaviors outside the 
interval —1.5 ^ s < 1.5; not so much for positive s~values (corresponding to the chiral 
expansion), but they differ quite dramatically for negative s-values (corresponding to the 

OPE expansion) where one finds: = 



while 



^WKiQ')U ..-^-^+'4 + o(^]. (4.13) 



2 ^^y^ ^IHM,_,^ 3^2 ^4 ' ^6 ' U8 



We can see that, except for the first term, the signs are different and the discrepancy in 
magnitude increases for each successive term. Yet, the shapes of the two correlation functions 
in the Euclidean, as seen in Fig. 5, is not so different. 

From these considerations we conclude that it is the Mellin transform of the spectral 
function (and its analytic continuation) which encodes all the details of the underlying dy- 
namics. One can also see that the determination of higher order condensates, the equivalent 
of the coefficients in the high-Q^ expansion in Eqs. (I4.12p and (I4.13P is very sensitive to the 
nature of the hadronic spectrum. This is an important issue for phenomenology because, for 
example, the coefficient of 0{1/Q^) in the expansion of nLR(Q^) plays a crucial role in our 
understanding of direct CP-violation in iT-decays, the so called e'/e contribution [30]. Un- 
fortunately, this 0{1/Q^) term, as well as the higher order ones, appear to be very sensitive 
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to the underlying dynamics. It is not surprising that the phenomenological determinations 
of these coefficients found in the hterature differ so much. 

V The Riemann Zeros and Sum Rules 

In the previous section we have considered a model from Analytic Number Theory to illustrate 
some physical features of Large-Nc QCD which we would like to understand. Here I shall 
do the contrary, I will consider the physical framework of Dispersion Relations in Quantum 
Field Theory to discuss a well known problem in Mathematics. The problem has to do with 
the positions of the zeros of the Riemann zeta function defined by Eq. 13.21 and its analytic 
continuation El which extends to all values of s except at s = 1 where it has a simple pole 
with residue 1. 

The interesting object for our purposes is the logarithmic derivative of the Riemann zeta 
function which, using the Euler product expression in Eq. 13.21 can be written as a Dirichlet 
Series: 

-7^= E iog(P)i;^-" 

primes p k=l 
oo 

= A(n) n-' , Re(s) > 1 , (5.1) 

n=l 

^ V ' 

Dirichlet Series 

where the A(n) (n integer) are the so called Von Mangoldt Amplitudes: 

^^^^ ^ , otherwise (^•^) 

A plot of the Von Mangoldt amplitudes for the first hundred integers is shown in Fig. 6. In 
^^See e.g. refs. 41, 42 and references therein. 

^^For a nice elementary treatment of the Riemann zeta function see e.g. refs. [431 144) . 
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20 40 60 80 100 

Fig. 6 The Von Mangoldt values A(n) for the first 100 integers. The upper curve corresponds to 
integer values coinciding with a prime number. The background points correspond to integers which 
are more than a power of a prime number. They appear as successive horizontal points. 

fact, there exists an explicit expression for the analytic continuation of the Dirichlet series in 
Eg. IS-H which follows from Hadamard's product formula for C(s) |44] : 

The poles at s = —2, —4 — 6 • • • correspond to the trivial zeros of C,{s) and the poles at s = p 
to all the remaining zeros. The non-trivial zeros satisfy < Re(p) < 1 , and because of the 
Symmetry Relation s — t- 1 — s, they must be located symmetrically relative to the vertical 
line Re(s) = 1/2, the so called critical line. The famous Riemann Hypothesis (RH) states 
that all the non-trivial zeros have Re(s) = 1/2 . Numerically, all the non-trivial zeros which 
have been evaluated so far do indeed satisfy the RH. 

V.l Question 

What are the properties of a Large-^c QCD-like Green's Function which has as a spectral 
function the one associated to the Von Mangoldt Dirichlet Series i.e., 

1 °° 

-lmnvonM(i) = y A(n) nM^ 5{t - nM^) , (5.4) 

n=l 

with the ampltudes A{n) given in Eq. ( 15. ^)) and plotted in Fig. 6 ? 

I wish to clarify the meaning of this question. From the point of view of QCD, this spectral 
function can only be considered, at best, as a toy model of large-Nc (perhaps as a toy model 
of duality violations). What seems interesting to me, however, is the fact that such an ab- 
stract mathematical question as the RH can be phrased, as we shall see next, in terms of the 
language familiar to physicists working in Quantum Field Theory. 
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Let us call HvonMl'?^) the two-point function which has as an imaginary part the Von 
Mangoldt spectral function in Eq. 15.41 The function nvonM('i'^) then obeys a dispersion 
relation modulo subtractions. The fact that the Van Mangoldt Dirichlet series in Eq. 15.11 is 
defined for Re(s) > 1 fixes the number of the required subtractions. They can be removed 
by taking two derivatives in the dispersion relation. This defines a function 'PvonM(-2) in the 
euclidean {z = '^), which is the analog of Eq. (j2.3p : 

'P\onu{z) = 7^7^2\2-'^VonM(<5 ) = 2 / dt — ImnvonM(i) 

[OLJ ) Jq (t + j it 

oo 

= 2^A(n) (5.5) 

n=l ^ ' 

with the corresponding Mellin-Barnes representation (c = Re(s) G ] + 1, +2[) 

c+«oo 

T'vonMl^) = j ds SvonM(s) T{s + l)r(2 - s) , (5.6) 

c—ioo 

where SvonM('5) is given in Eq. 15.31 Everything is explicitly known in this representation 
and we can now apply the Inverse Mapping Theorem of ref. [32J to compute the asymptotic 
behaviour of PvonM(-z)- 

The interesting expansion is the short distance expansion corresponding to large 
(small-z). This expansion is governed by the singularities at the left of the fundamental 
strip i.e. s < 1. The singularities at s = —1,-2,— 3, ••• generated by the r(s + 1) factor 
in the integrand and by the poles in SvonM('S) at s = —2, —4, —6, • • • corresponding to the 
trivial zeros of C(s) give rise to odd powers of 

. 2X 2n+l 

' n = 0,l,2,3--- , (5.7) 



as well as to even powers (n = 1, 2, 3 • • • ) of 

/M2\^" 02 /m2\^" 

These terms are rather analogous to the usual power terms which originate in the OPE in 
Quantum Field Theory. The leading singularity at s = 1 gives rise to the leading asymptotic 
behaviour for large, which is: 

^ - r(2)r(i).-i = ^ , (5.9) 



s-l ^ ^ ^ ^ M"^ 

in fact rather similar to the leading behaviour of a QCD-like two-point function generated 
by a scalar current. 

The interesting terms are of course the ones generated by the next -to-leading singularities 
at s = p in the SvonMls) function; i.e. the ones induced by the non-trivial zeros of the 
Riemann zeta function. They give rise to non-power terms: 

-^^-^(p + l)^(2-p)z-^ (5.10) 
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which appear in pairs of O (^%^^ ^ (^) ' ^o^^^^^^*^ ^ly an oscihating behaviour 

in Q^. In the particular case where p = 1/2 ±17], i.e. for the zeros satisfying the RH, these 

terms collapse to a unique non-power behaviour of O 

-1 -1 
+ 



(i + ir]) s - (5 - iv) 

0^1+ V vr / 

cos Ulog-Tg- , (5.11) 



M2 2 coshvTT? V Q"^ 

modulated by oscillating cos ^r/log factors, one for each value of r] along the critical line 
of zeros, with amplitudes ^~^2^ coshTrr? '^hich decay exponentially for rj large. 

Conclusion: The Riemann Hypothesis is equivalent to the existence of a unique type of 
non-power terms of O ^ in the short distance expansion of the two-point function 

associated to the Von Mangoldt spectral function in Eq. \ 5.4\ 

V.2 Quantum Mechanics Sum Rules 

It is perhaps helpful to discuss the previous considerations within a more general Statistical 
Mechanics Framework. The discussion applies to Hamiltonians Ti with no explicit time 
dependence. 

The probability transition amplitude in Quantum Mechanics is defined as 

{qf,tf\quti) = {qf\e''^^'f-'^'>\qi). (5.12) 

The evolution in imaginary time leads to a Statistical Mechanics interpretation which is 
characterized by the Partition Function 

Z = Tvex.p-l3'H. (5.13) 

With tf — ti = —i(3, we then have the spectral representation: 

{qf\e-'nq^) = V e^^^" {qf\n){n\q,) . (5.14) 




In particular, the lim/3 — c«; i.e., T — 0, [/3 = is governed by the ground state 

contribution: 

lim {qf\e-^^) ~ e-^^o^-o (9^)^2(9/) • (5-15) 

/3— !>oo 

In general: 

Eo= lim -^logTre"'^^. (5.16) 

/3->-+oo p 

The following quantity 

M{^) = {qf = 0\e-^^\qi = 0) = |^„(0)|2 exp[-/?S„] . (5.17) 
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is of special interest to us because it provides the Quantum Mechanics framework for the 
equivalent of the Quantum Field Theory Sum Rules. 

The relevant Hamiltonian for our purposes is one which has levels 



En = nEo 



(5.18) 



and wave functions at the origin 



|^'„(0)|2 = A(n), 



(5.19) 



i.e. the Von Mangoldt amplitudes defined in Eq. ()5.2p . Inserting this ansatz in Eq. (|5.17p . 
the corresponding Mellin-Barnes representation is then 



c+ioo 

-MvonM(/3) = ^ / dsi^Eo)-' SvonM(s) r(s) , 



(5.20) 



with a fundamental strip defined now by the interval c = Re(s) G ] + 1, oo[ and SvonM(s) the 
same expression as in Eq. 15. 3[ The inverse mapping theorem applied to this formula gives us 
the expansion at small /3; i.e. the expansion at high temperature (/3 = 



M 



VanM 



(/3) 



1 



PEo 



= V [r(l/2 + ir?)e*^''^°s(/3So) + r(l/2 - ir^)e-''i^°s{l3Eo) 
3En — ^ L 



+ logi^ + /3i.o(log(8vr)-l + ^-^) 



+ {f3Eof 
+ 0[{(3Eo 



log(/3i^o)-^-^log(4vr) + i7E-E;^ 



(5.21) 



In this expression the leading behaviour in the first line is the one associated to the trivial 
singularity at s = 1, while the second term in the first line gives the asymptotic behaviour 
reflected by the non-trivial zeros located at s = 1/2. The Riemann Hypothesis implies that 
these are the only possible non-power terms in the expansion. A constant term as well as 
odd power terms in {(3Eq), with the leading one shown in the second line, are generated 
because of the T{s) factor in Eq. 15.201 with coefficients which, except for the constant term, 
depend on the location of the non-trivial zeros (the sum over p's). The trivial singularities at 
s = —2n generate powers in (/3i?o)^" (the leading one is shown in the third line) modulated 
by a log(/3So) factor and powers of (PEq)'^^. The latter are modulated by coefficients which 
also depend on the positions of the non-trivial zeros (the sum over /)'s). 

An interesting question which I have been investigating in collaboration with Josep Taron 
is: what is the equivalent potential which produces a spectrum of discrete levels En = tiEq 
with corresponding wave functions at the origin known in modulus |^'n(0)p = A(n) (i.e. the 
Von Mangoldt amplitudes defined in Eq. i5.^} )? I hope to have an answer for Raymond's 
next birthday... 
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